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Survival and extinction results for a patch 
model with sexual reproduction 

Eric Foxair and Nicolas Lanchier^ 

Abstract This article is concerned with a version of the contact process with sexual repro¬ 
duction on a graph with two levels of interactions modeling metapopulations. The population 
is spatially distributed into patches and offspring are produced in each patch at a rate pro¬ 
portional to the number of pairs of individuals in the patch (sexual reproduction) rather than 
simply the number of individuals as in the basic contact process. Offspring produced at a 
given patch either stay in their parents’ patch or are sent to a nearby patch with some fixed 
probabilities. As the patch size tends to infinity, we identify a mean-field limit consisting of an 
infinite set of coupled differential equations. For the mean-field equations, we find explicit con¬ 
ditions for survival and extinction that we call expansion and retreat. Using duality techniques 
to compare the stochastic model to its mean-field limit, we find that expansion and retreat are 
also precisely the conditions needed to ensure survival and extinction of the stochastic model 
when the patch size is large. In addition, we study the dependence of survival on the dispersal 
range. We find that, with probability close to one and for a certain set of parameters, the 
metapopulation survives in the presence of nearest neighbor interactions while it dies out in 
the presence of long range interactions, suggesting that the best strategy for the population 
to spread in space is to use intermediate dispersal ranges. 


1. Introduction 

The term Allee effect refers to a certain process that leads to decreasing net population growth 
with decreasing density [T]. In case the growth rate becomes negative at low density, this monotone 
relationship results in the existence of a so-called Allee threshold below which populations are at 
high risk of being driven toward extinction. This phenomenon may be due to various ecological 
factors: failure to locate mates, inbreeding depression, failure to satiate predators, lack of coopera¬ 
tive feeding, etc. Research on this topic is copious and is reviewed in [3| but rigorous mathematical 
analyses of stochastic spatial models that include an Allee effect are much more limited. 

In the model proposed by Borrello [3], each site of the infinite regular lattice represents a patch 
that can host a local population, and a strong Allee effect is included in the form of a varying 
individual death rate taking a larger value for local populations below some threshold. The model 
is used to show that, when only small flocks of individuals can migrate from patch to patch, the 
metapopulation goes extinct whereas survival is possible if large enough flocks of individuals can 
migrate. The framework of interacting particle systems has also been used in BM to study the 
consequence of an Allee effect. There, the model is a modification of the averaging process that also 
includes a threshold: local populations below this threshold go extinct whereas local populations 
above this threshold expand to their carrying capacity, each at rate one. The key component in 
this work is the topology of the network of interactions rather than the size of the migrating flocks, 
and the analysis of the process starting from a single occupied patch on various graphs indicates 
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that the probability of long-term survival of the metapopulation decreases to zero as the degree 
of the network of interactions increases to infinity. This result suggests that long range dispersal 
promotes extinction of metapopulations subject to a strong Allee effect. 

The modeling approach of the present paper is somewhat different. The model we propose is a 
version of the contact process with sexual reproduction mm on a graph that includes two levels 
of interactions modeling metapopulations: individuals are produced within each patch at a rate 
proportional to the number of pairs of individuals in the patch rather than a rate proportional to 
the number of individuals in the patch. This birth mechanism (sexual reproduction) reflects the 
difficulty to locate mates in patches at low density, which has been identified by ecologists as one of 
the most common causes of Allee effect. In particular, while the Allee effect is forced into the model 
in the form of a threshold parameter in lasi, it is on the contrary naturally induced by the birth 
mechanism in the model considered in this paper. A useful consequence of this fact is both our model 
and its mean-field equations have a dual process, moreover the stochastic model’s dual process is 
a straightforward truncation of the mean-field’s dual process. With this relationship in hand, we 
can easily estimate the occupation density of each patch in the stochastic process by the solutions 
to the mean-field equations, and using block constructions, we can show the following. When a 
localized population spreads in the mean-field equations, we have survival of the stochastic process 
and existence of a non-trivial stationary distribution. When a localized low-density region erodes 
the surrounding population in the mean-field equations, the same occurs in the stochastic process 
and implies weak convergence to the vacant configuration. We analyze survival and extinction 
starting from either a large block of occupied patches, or a single occupied patch, and also assess 
the dependence of survival on the dispersal range. 

2. Model description and main results 

Each integer x G Z represents a patch that can host a local population of up to N individuals, 
and we think of each site as a set of spatial locations that can be either empty or occupied by one 
individual. By convention, we use bold letters to denote these spatial locations: 

X := (x,j) e Dn :=Zx {1,2,...,N} 

are the possible spatial locations at site/patch x. The model we consider is a continuous-time 
Markov chain whose state at time t is a spatial configuration 

rjt : Dj\f —{0,1} where 0 = empty and 1 = occupied. 

Each individual dies at rate one. Offspring produced at a given patch are sent to a spatial location 
chosen uniformly at random from either the parents’ patch or a neighboring patch. In either case, 
offspring are produced at a rate proportional to N times the fraction of pairs of spatial locations 
which are occupied and the birth is suppressed if the target is already occupied. The proportionality 
constant is denoted by a for offspring sent within their parents’ patch and by b for offspring sent 
outside their parents’ patch. Sexual reproduction is modeled by the fact that the birth rate is related 
to the number of occupied pairs and the reason for multiplying by N is to have births and deaths 
occurring at the same time scale. For simplicity, we consider symmetric interval neighborhoods 
but most of the basic facts, and some of the estimates, are true for any translation-invariant 
neighborhood, on Z or on more general transitive graphs. More formally, for x, y G Z we write 

X ~ y if and only if x ^ y and \x — y\ < M 
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where M is the dispersal range, and dehne the projection map 


vr : X := {x,j) G Djy i-)- 7r(x) := x G Z. 

For all X G -DAT, we also dehne the sets 

A(x) := set of potential parents’ pairs within the patch containing x 
:= {(y, z) G -Dtv X F>Ar : y / z and 7r(x) = 7r(y) = 7r(z)} 

B{x.) := set of potential parents’ pairs near the patch containing x 

:= {(y, z) e Dn X Dn ■■ y ^ z and 7r(x) ~ 7r(y) = 7r(z)}. 


The dynamics is then described by the Markov generator 


L-fiv) = [/(^x,o) -/(h)] 

X 


+ E 


\2M N{N- 


(y.z)GA(x) 


( 1 ) 


jy h(y)h(z)) [/(hx,i)-/(h)] 


(y,z)eS(x) 


where conhguration is obtained from 77 by setting the state at x equal to i. We shall call this 
process the microscopic representation. To study this model, it is convenient to also consider its 
mesoscopic representation that keeps track of the metapopulation at the patch level rather than at 
the individual level. This new process is simply obtained by setting 


6 (a:) := Ex:7r(x)=x hf(x) for all {x,t) G Z x M+. 

In words, the process counts the number of individuals at each patch. Note that, since the particular 
locations of the individuals within each patch is unimportant from a dynamical point of view, this 
new process is again a Markov process, and its Markov generator is given by 


L+fiC) = 
+ 

+ 


E ^(^) [/(^-) - /(^)] 

X 

E n{N- 1) - 1)(A - ((x)) l/(U+) - /({)! 

E E NiN- 1) [/&+)-/({)] 

X y'^X ^ ' 


( 2 ) 


where conhguration is obtained from ^ by adding/removing one individual at x. The analog 
of this model derived from the basic contact process rather than the contact process with sexual 
reproduction has been studied in [ 2 ] where it is proved that 


• the process survives when a + 6 > 1 and N is sufficiently large, 

• the process dies out for all values of the parameter N when a + h < 1. 


The analysis of the stochastic process ©-(H is more challenging due to the complexity of the birth 
mechanism. Our approach is to use duality techniques to show that, at least in bounded space-time 
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regions and when N is large, the stochastic process can be well approximated by the system of 
coupled differential equations for u = {ux)x&-, called the mean-field equations, given by 

K = (aul + Ux) - Ux for all x € Z. (3) 

y'^x ^ 

The long-term behavior of the process can then be deduced from properties of the mean-field 
equations combined with block constructions. 

The mean-field equations - Starting from u(0) constant, the profile u{t) remains constant across 
space at all times and solves the single differential equation 

u' = ru^ (1 — u) — u where r := a -\-b. (4) 

Some basic algebra shows that 

• when r < 4, there is a unique equilibrium, namely 0, 

• when r = 4, we have the pair of equilibria: 0 and 1/2, 

• when r > 4, there are three equilibria: 0, 

U-:=l/2 — w and u+:=l/2-\-w where tc = (1/4 — 1/r)^/^ 
with 0 and u+ stable and U- unstable. 

This motivates the following definition, which is the key to understanding the system ([3|) starting 
from more general profiles: for r = a -|- 6 > 4, we say that 

• expansion occurs if there is u with tt_ < tt < so that 

u(0) = u l{x < 0) implies that ui(to) = u for some to > 0; 

• retreat occurs if there are and u* with 0 < u* < ri_ < < ti* so that 

u(0) = u* l{x < 0) -|- tt* l(x > 0) implies that U-fitfi) = u* for some to > 0. 

Notice that retreat is also defined for r = 4 by letting rt_ = = 1/2. Moreover, if r < 4 we shall 

say that retreat occurs, since then the conclusion of the upcoming Theorem [1] holds for any > 0. 
Expansion implies that starting from a large enough occupied block of patches, the population 
spreads at a linear rate. Retreat implies the opposite: starting from a large enough vacant block of 
patches, the population dies out within a region that grows linearly in time. This is summarized by 
the following result, which can be deduced from more general results of m- We give a more direct 
proof in Section [3 using our duality theory. 

Theorem 1 - We have the following implications for expansion and retreat. 

• If expansion occurs, there are u, L, xq, 5, c > 0 with u- < u < so that 

i^a:(0) > u for all |x| < L implies that Ux{t) > u-\- 5 for all \x\ < ct — xq. 

• If retreat occurs, there are u, L, xq, 5, c > 0 with 0 < u < U- so that 

^^a;(0) < u for all |x| < L implies that Ux{t) < u — 6 for all \x\ < ct — xq. 
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It is not possible to obtain exact threshold values for expansion or retreat. However, we can obtain 
some decent estimates, which among other things implies that when M = 1, retreat occurs for an 
open set of values of (a, 5) satisfying a + 6 > 4. 

Theorem 2 - Expansion and retreat are open conditions in that the sets 

{(a, 6) : expansion occurs} and {{a,b) : retreat occurs} 
are open subsets of x Moreover, if M = 1, then 

• expansion occurs when a + 6/2 > 4 and b > 8/9, 

• retreat occurs when a + 6 < 4 and 6 > 0. 

Section [7] also collects results for survival of the system of ordinary differential equations starting 
with a single fully occupied patch and all the other patches empty. 

The stochastic process - We now state our results for the stochastic process. Constructing 
the system graphically from a collection of independent Poisson processes and using standard 
coupling arguments, one easily proves monotonicity with respect to the birth rates a and 6, as well 
as attractiveness. Note however that basic coupling arguments do not imply monotonicity with 
respect to the patch capacity N or the dispersal range M. Attractiveness implies in particular 
that the limiting distribution of the process starting from the all occupied configuration exists, so 
survival and extinction can be studied through this invariant measure looking at whether it is a 
nontrivial distribution or the point mass at the all empty configuration, that we denote from now 
on by 0. Using duality techniques and block constructions as well as properties of the system of 
differential equations given by Theorem [2l we can prove the following two results. 

Theorem 3 - Assume that expansion occurs in the mean-held equations. Then, for all N large, 
the stochastic process has a nontrivial stationary distribution, and there is L such that 

limjv^cxD P = 0 for some t > 0 | ^o{x) = N for all x G [—L, L]) = 0. 

Theorem 4 - Assume that retreat occurs in the mean-held equations. Then, for all N large, the 
stochastic process converges in distribution to the point mass at 0. 

To answer an important ecological question, namely, whether an alien species established in one 
patch either successfully spreads in space or is doomed to extinction, we now study the probability 
of long-term survival for the process starting with a single fully occupied patch, that is the process 
starting from the initial configuration 

^o(O) = N and ^ 0 ( 2 ;) = 0 for all x yf 0. 

From now on, we let Pq denote the law of the process starting from this conhguration. Relying 
again on properties of the system ([3]) together with duality techniques and block constructions, we 
obtain the following sufficient condition for successful invasion. Note the conditions on the birth 
rates are slightly stronger than for expansion, since we want the population to spread starting from 
a single occupied patch and not just from a sufficiently large finite block. 

Theorem 5 - Assume that M = 1. Then, 

limAr_).oo Po {Ct = 0 for some t > 0) = 0 whenever a -\- 6/2 > 4 and 6 > 2. 
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Finally, we study he dependence of survival on the dispersal range M. Theorem [4] shows that 
extinction occurs starting from any configuration when the birth parameters are small enough and 
the patch capacity is large. In contrast, our last theorem focuses on the process starting with a 
single fully occupied patch, in which case, regardless of the birth rates and the capacity of the 
patches, extinction occurs with probability close to one when the dispersal range is sufficiently 
large. 

Theorem 6 - For all M large, we have 

Po {Ct (1 for all t > 0) < {1/2 + bN{l — a/4)~^) when o<4 

< (1/2 + (6/2)(A^ + 2)^) when a = 4 

< (1/2 + 6 (a/4 — 1)“^ (a/4)'^+^) when a > 4. 

There are three different estimates because the survival probability is related to the time to extinc¬ 
tion of a patch in isolation, which scales differently depending on whether the inner birth rate is 
subcritical, critical or supercritical. In either case, the theorem shows that the survival probability 
decreases to zero as the dispersal range increases to infinity, indicating that long range dispersal 
promotes extinction of metapopulations subject to a strong Allee effect caused by sexual reproduc¬ 
tion. In particular, the effects of dispersal are somewhat opposite for the process with and without 
sexual reproduction since, as proved in [2], in the presence of long range dispersal, the basic contact 
process approaches a branching process with critical values for survival significantly smaller than 
that of the process with nearest neighbor interactions. 

The paper is laid out as follows. In Section [3l we introduce the dual to the stochastic model and 
show monotonicity of both. In Section [H we introduce the dual to the mean-field equations and 
use it to prove some useful approximation results for solutions to the mean-field equations, that 
can be viewed as localized and quantitative statements of continuity with respect to initial and 
boundary data. In Section [5l we show how to obtain the dual to the stochastic model by truncation 
of the dual to the mean-field equations, and conclude the two duals coincide so long as there is no 
collision of particles. In Section [H we use the above agreement of duals to compare the occupation 
density at each patch in the stochastic model to the mean-field value, with an error term that 
is proportional to the collision probability. In Section [7l we pause for a moment to establish the 
stated properties of the mean-field equations. This uses only the results of Section [H In Section [HI 
we prove Theorem [3l using a block construction and the occupation density estimates of Section El 
In Section [9l we prove Theorem 01 using the same idea, but with larger blocks and with more care 
to ensure the establishment of a completely vacant zone that grows over time. In Section m we 
prove Theorem El in the same manner as Theorem El Finally, in Section fTTl we show that, with high 
probability when M is large, the initially occupied patch dies out before any two individuals born 
at that patch are sent to the same neighboring patch in order to obtain Theorem El 

3. The dual process of the patch model 

To prove Theorems EHSl we need to approximate the stochastic process using the mean-field equa¬ 
tions, and we achieve this by constructing dual processes for both systems, which we then relate 
to one another. We refer to the stochastic process with patch size N as the N-patch model, and 
to its dual process as the N-dual. This will help to distinguish these dual processes from the dual 
of the mean-field equations introduced in the next section. The main objective of this section is 
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to construct both the A^-patch model and its dual process from a graphical representation. Then, 
using this graphical representation, we prove that both processes are monotone with respect to the 
inner and outer birth rates a and b. To define the graphical representation, for each microscopic 
spatial location x, we introduce the following random variables. 

• For each (y,z) € ^(x), 

{an(x, y,z) : n > 0} := Poisson point process with rate a{N{N — 1))“^. 

• For each (y,z) € B{x), 

{6„(x, y, z) : n > 0} := Poisson point process with rate (6/2M) {N{N — 1))“^. 

• In addition, {(in(x) : n > 0} := Poisson point process with rate 1. 

The microscopic process is constructed from these Poisson processes as follows: 

• Births: at time t = an(x,y, z) or t = 6„(x, y,z), we set 

774 (x) = 1 when 774 . (y) 774 . (z) = 1 

= 774 - (x) when 774 _(y) r 74 -(z) = 0 . 

• Deaths: at time t = dn{^), we set 774 (x) = 0. 

We also construct fis{'w,t), the dual process starting at space-time point (w,t), to keep track of 
the state at this point based on the configuration at time t — s. The dual process of the contact 
process with sexual reproduction consists of a collection of hnite subsets of D^, the set of spatial 
locations. To describe its dynamics, it is convenient to introduce 

:= {x € Dw : X € B for some B € fjs{w,t)} = UiJGiiAw.t) ^ 

which we call the active set. 

1. The process starts from the singleton 77 o(w, t) = {{w}}. 

2. Births: if site x G A^_ where either 

s = t — an(x, y, z) or s = t — 6n(x, y, z) for some n>0 

then, for each set B G 77 s_(w,t) that contains x, we add the set which is obtained from B by 
removing x and adding its parents’ sites y and z, i.e., 

%(w,t) := fis-(w,t) U {(B - {x}) U {y,z} : X G B G 77s_(w,t)}. 

3. Deaths: if site x G A^_ where s = t — d„(x) for some tt, > 0 then we remove from the dual 
process all the sets that contain x, i.e., 

77s(w,t) := 77s_(w,t) - {B G r7<i-(w,t) : X G B}. 

See Figure [U for a picture. The dual process allows us to deduce the state of site w at time t from 
the configuration at earlier times. More precisely, identifying the microscopic process with the set 
of occupied sites, the construction of the dual process implies the duality relationship 

w G 774 if and only if B C r/ts for some B G 77s(w, t). (5) 

Using a basic coupling argument, we now prove that the Wpatch model and the A^-dual are 
monotone with respect to the inner and outer birth rates. 
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(x, t) 



patch X — 1 patch x patch x + 1 patch x + 2 


Figure 1. Picture of the A'^-dual for N = 5 with time going up. Arrows represent birth events with the offspring at 
the head of the arrow and the parents’ location being the two dots at the tail of the arrow, while crosses represent 
death events. The bold lines refer to the active set. In our realization, the dual at the bottom of the picture consists 
of only one set of cardinal five, so site x is occupied at time t if and only if these five sites are occupied. 


Lemma 7 - The processes rjt and ris(w,t) are nondecreasing with respect to a and b. 

Proof. Fix ai < a 2 and 6 i < 62 and, for i = 1,2, let 

rfi := the A^-patch model with inner and outer birth rates ai and hi. 

Then, construct the first A^-patch model rjl from the graphical representation above with a = oi 
and 6 = 61 . Basic properties of Poisson processes imply that the process can be constructed from 
the same graphical representation supplemented with additional independent Poisson processes 
with intensity 02 — oi for inner births and 62 ~ for outer births. This defines a coupling of the 
two A^-patch models for which one easily shows that 

P {rjj C r]f) = 1 for all t > 0 whenever C rjQ. 

This shows that the Af-patch model is nondecreasing with respect to the inner and outer birth 
rates. The monotonicity of the Ai-dual can be proved similarly. □ 

4. The dual process of the mean-field equations 

We now discuss the dual of the mean-field equations, which we call the limiting dual because it 
appears to be the the limit in distribution of the Af-dual, and denote by Ct- Its state space consists 
of the set of finite collections of finite sets of points in := Z x (0,1). As previously, the dynamics 
is described using the active set given by 

At := {(x, w) (i D : {x,w) (i B for some B G Ct} = Useft 

Then, the process Q has the following transitions. 
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• For each {x, w) G At, at rate a, add to Q all sets obtained from sets B G Q such that {x, w) G B 
by removing {x, w) and adding {x, Wi), (x, W 2 ) where Wi,W 2 are independent uniform random 
variables on the interval (0,1). 

• For each {x,w) G At and each y ~ x, at rate b/{2M), add to C,t all sets obtained from 
sets B G Ct such that (x,u;) G B hy removing (x,lc) and adding {y,wi), {y,W 2 ) where wi,W 2 
are independent uniform random random variables on (0,1). 

• For each {x,w) G At, at rate 1, remove from Q all sets B containing {x,w). 

We now exhibit the connection between the limiting dual and the mean-field equations. Suppose u 
is given. For B G (t, writing B = {(xijtui),..., {xk,Wk)}, say that 

B is good for u if and only if Wj < Uxj for all j = 1, 2,... , A: 

and write Ct ~ u if there exists B G Ct that is good for u. Letting w ~ Uniform (0,1), for t > 0 

and u G Z9, define the function cj) : M_|_ x /C —)• /C where /C = [0,1]^ by 

{M^))x = PiCt ^Co = {{ix,w)}}). (6) 

We want to show that cjit gives us the solutions to the mean-field equations. Our first task is to 
show that it has the semigroup property. 

Lemma 8 - For (p as defined in ([6]) and s,t > Q, we have (pt+s = <Pt° 4’s- 

Proof. For an illustration of what is introduced in the proof, we refer to Figure [2j First, we 

re-express Ct as a labelled set of points in D, namely as {It, it), where 

• for each t > 0, the influence set It is a finite set of points in D and 

• for each t > 0, the labelling it : It ^ keeps track of the child-sibling-parent relation 
whenever the influence set It branches. 

The influence set starts at Iq = {(x,r(;)} and has the following transitions. 

• For each {x,w) G It, at rate a, add to It the pair of points (x,ri;i) and {x,W 2 ) where wi,W 2 
are independent Uniform (0,1) random variables. 

• For each (x, w) G It and each y ~ x, at rate b/ (2M), add to It the pair of points (y, wi), (y, W 2 ) 
where wi,W 2 are independent Uniform (0,1) random variables. 

• For each {x,w) G It, at rate 1, remove {x,w) from fl. 

Since newly added points {x,w) have w ~ Uniform (0,1), with probability one they do not collide 
with existing points, so is a branching random walk. Also, we have At C It under the obvious 
coupling, and strict inclusion is possible since the removal of a point from the influence set causes 
the removal of all sets B G Ct that contain this point. 

The labelling is then defined as follows. Each point in gets three labels: its own personal 
label, a generation number, and the label of its parent, with the parent label of the initial point 
being unimportant. If a point is removed from It, its personal label and generation number are 
not reused, and the labelling of a given point does not change over time. Personal labels and 
generation numbers are in increasing order of appearance, so the single point in Iq gets the label 1 
and generation number 1, and if personal labels !,...,& and generation numbers 1,..., m are in 
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( 1 , 1 ) 



Figure 2. Picture of the limiting dual with time going up. Arrows, dots and crosses have the same interpretation as 
in Figure [T] while the bold lines now refer to the influence set. The pairs of numbers in the picture are respectively 
the personal label and generation number of each of the sites that are included in the influence set. 


use and a new pair of points appears, then they are given personal labels k + 1 and A:+ 2, generation 
number m + 1 , and the label of their parent. 

Given {Is,£s)o<s<t, we now identify a set of labels active at time t recursively as follows. For 
each (x, w) € It, if tu < Ux, the label of (x, w) is active. If there are two active labels with the same 
generation number, their common parent label is active. Since the generation number of a parent is 
smaller than that of its child, and since, with probability one, after a hnite time only a finite set of 
generation numbers have been assigned, if we repeat the last step, after a finite number of iterations 
all active labels have been found. The intuition behind this construction can be understood returning 
to the microscopic process as follows: whenever there is a birth event {y, z} —)• x, the target location 
becomes occupied if both parent’ locations are occupied. It follows by inspection that 

~ u if and only if label 1 is active at time t. 

We now observe the labelled influence set has the following composition property. Given s,t>0 
and {{Ir,ir) ■ 0 < r < t}, each particle in It evolves independently and in the same way as a single 
particle at time 0 Therefore, we can construct {{It+r,^t+r) : 0 < r < s} by first conditioning on It 
and then appending, to each (x,w) G It, an independent copy of the process started from {x,w). 
To maintain consistency with It it then suffices to re-label each copy while preserving the relation 
between the labels. The set of values {w : {x,w) G It} are independent Uniform (0,1) random 
variables, so it follows that the labels of points (x, w) G It are active at time t + s independently 
with probability Ux{s) = {4>s{u))x. Using the same fact, this can be equivalently stated by saying 
that the label of {x,w) G It is active if ru < Ux{s). Using the definition of 4>, we find that the label 
of the single point {x,w) G Iq is active with probability {(I)ti4‘s{'^)))x- Since it is also active with 
probability ( 04 +<j(u)) 3 ;, this completes the proof. □ 
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We say that (pt is the flow corresponding to Q if for each u € /C, the function u(t) = (/>t(u) 
is the unique solution to ([3]) with u(0) = u. 

Theorem 9 - The function (pt is the how corresponding to ([3]) . 

Proof. We need to check that 

dt(pt{rT) = F{u) for all t > 0 and u G /C (7) 

where F{u) is the right-hand side of Q. In light of Lemma [51 it suffices to show that equation ([7]) 
holds when t = 0 and u G /C, since 

dt(pt{^) = iimft^o+ - </>t(u)] 

= limh^Q+ h-^ [4>h{(pt{^)) - 

and (pt{rL) G fC. In order to prove dTj), we first observe that, given Co = and h> 0 small, 

the list of all possible events for Cft along with their probabilities are 

• El := {C/i = {{(x,!(;)}, {(x, rci), (x, 102 )}}} for some independent uniform tui, t(; 2 , which occurs 
with probability ah o{h), 

• £^ 2 ,y := {C/i = {{(^) ii')}! {( 2/5 ii'i)! ( 2 /) ii' 2 )}}} for each y ~ x and for some independent uni¬ 
form wi,W 2 , which occurs with probability bh/{2M) + o{h), 

• E 3 := {C/i = 0 }, which occurs with probability h -|- o{h), and 

• £'4 := {C/i = Co}) which occurs with probability 1 — (1 -|- a -|- b)h -|- o{h). 

In addition, since w,'Wi,W 2 are Uniform (0,1), 

P (C/I ~ u I £ 1 ) = P (to < Ux) + P{w > Ux, Wl < Ux, W 2 < Ux) = Ux + (1 - Ux) ul, 

P (C/I ~ u I E 2 ,y) = P {w < Ux) + P (w > Ux,Wi < Uy,W 2 < Uy) = Ux + {1 - Ux) Uy, 

-P (C/i ~ U I £ 3 ) = 0, 

P (C/i ~ U I £ 4 ) = P (to < Ux) = Ux, 

P (Co ~ u) = P (rc < Ux) = Ux. 

Putting things together and noting the cancellation of Ux [1 — (a + b)h], 

<9^0^ (u) I= lim;,^o+ [P iCh ~ u) - P (Co ~ u)] 

= a P (C„ ~ u I £i) + 6 (2M)-i P iChn \ £ 2 ,^) 

- (I + a + 6 ) P (C/i ~ u I £ 4 ) 

= a {I - Ux) ul + b {2M)-^ Y)yr^x(^ - '^x) y-l - Ux = £(u) 
from which the desired result follows. □ 

Theorem [9] implies some useful properties for the mean-field equations. The following is a direct 
consequence of Theorem [9] and the definition of (pt in terms of Ct ■ 

Corollary 10 ~ Let u{t) and v(t) be solutions to ([3]) with respective parameter values (au,bu) 
and (uyjby). Then, if a„ < a„ and by < by, 

u(0) < v(0) implies that u(t) < v(t) for all t > 0. 
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In the next lemma, we control the size of the influence set. 


Lemma 11 - Let := {x G Z : (x, w) € Ig} and 7 > 0. Then, 

P {Js C (—cxD, ct\ I Iq = {(0, Wo)}) > 1 — for all s < t and some c > 0. 

Proof. For 0 G M and t > 0, dehne 

mt( 0 ) := E I ^0 = {( 0 , wq)}) 

which encodes information about the distribution on Z of the points in 1^. By first conditioning on 
the value of It and noting that individual points evolve independently, we find 

dtmt{9) = mt{0)dsms{9)\^^^^ 

and we readily compute 

m := dsmg{9)\^^^ = a + (2M)-i " 1 

which, since mo{9) = 1, implies mt{9) = To control the spread of It, let 

Stic) ■■= \{{y,v) e Is : y > ct}\ = Yt{y,v)els'^iy > c G M. 

Since Stic) is integer valued and l(y > ct) < for each 9 G M+, 

PiStic)>0\ Jo = {( 0 , 1 x 0 )}) 

= EnGN* PiStic) = n\Io = {( 0 , 1 x 0 )}) < EiStic) I lo = {( 0 , 1 x 0 )}) ( 8 ) 

< ^(E(,,.)e/.e^"-'‘=M^o = {(0,ixo)}) = e-^'=‘m,(0). 

Now, we observe that, when b > 0, there exists 9 large such that ii9) > 0, in which case ms(0) is 
non-decreasing in s. In particular, using ([ 8 ]), we get, for c > ii9), 

P(5{(c) >0 I /o = {( 0 ,ixo)}) < e-^^*m,( 0 ) 

< mti9) < _ g-7(0)i 

for all s < t, where 7 ( 0 ) = 9c — ii9) > 0. On the other hand, when 6 = 0, the projection of the 
influence set on Z reduces to a singleton therefore 5{(c) = 0 and ([9]) is trivial. Finally, we use 
that 7 ( 0 ) ^ 00 as c —>■ 00 together with Q to conclude that, for all 7 > 0 , 

PiJs C (- 00 , ct] I lo = {(0,ixo)}) 

= P iStic) = 0 I Jo = {(0, ixo)}) > 1 - e-P for all s<t 


for some c > 0. This completes the proof. □ 

Lemma [11] has the following pleasant consequence for ([3|). 


Survival and extinction results for a patch model with sexual reproduction 
Lemma 12 - There exists c ,7 > 0 so that if u(t), v(i) are solutions of ([3]), 

'Wx(O) = i'a;(0) for all X G [y — ct,y + ct] 

implies that \uy{s) — Uy(s)| < 2 e“'^* for all s < t. 

In particular, for each €,t > 0, there is L > 0 such that 

Uxio) = Vxio) for all x G [y - L,y + L] 

implies that |rtj/(s) — '(;y(s)| < e for all s < t. 

Proof. Let A = {x G Tj : Ux = Vx} and s < t. Then, by definition of 4>s, 

P (0 ~ u I J, c = P (G ~ V I J, c A). 

Using Lemma [H] and (jlOh . we deduce that there exist c ,7 > 0 such that 

\uy{s) - Vy{s)\ = - ((?!)Gv))y| < P{JsnA'^^0\Io = {iy,w)}) 

< P{Js {-oo,y + cf\^ 0 \ Io = {{y,w)}) 

+ P (Js ^ [y - ct, 00 ) ^ 0 \ Io = {{y, w)}) 

< 2 P(Js G (-oojct] I Jo = {( 0 ,u;o)}) < 2 e“P 
which proves the first part. The second statement is an easy consequence of the first. □ 

In order to help with block constructions later on, we conclude this section with two more es¬ 
timates about truncated versions of our models. 

Definition 13 (restriction) - For integer K > 0 and fixed hut arbitrary boundary values Ux{t), 
K < \x\ < K + M, t > 0 that are continuous functions of t, we define the restrictions of the limiting 
dual and mean-held model ^Ux{t) as follows. 

• Given Co such that Jq C [—K,K], the restriction is obtained from G by freezing the 
evolution of any point {y, w) G It that lands in \y\ > K and declaring it good if w < Uy{t). 

• Given n(0), the restriction ^Ux(t) is dehned to be the solution to ([3|) where only the val¬ 
ues Ux{t) with Gl < k are determined by ([3|), and using the given boundary values Uy{t) 
if y ^ X and \y\ > K. 

The two most obvious choices of boundary values are Uy{t) = 0 and Uy{t) = 1, that we call the 
lower and upper boundary values. Dehning as before but in terms of h follows that 

'^0Gu(O)) = ^Ux{t) for all xG[—K,K] and t > 0. 

One verifies also that with lower boundary values, ^Co C Co implies ^Ct U Ct t > 0, and the 
same holds for upper boundary values but with the inclusion reversed. Since the event C ~ u is 
increasing with C in the sense that if C ~ u and C P then (^' ~ u, for lower boundary values we 
obtain that 

^u(O) < u(0) on [—K, K] implies ^u(t) < u{t) on [—K, K] for all t > 0 

and the same holds for upper boundary values but with the inequality reversed. The following result 
fills in the other side of the inequality in both situations. 
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Lemma 14 - Setting ^Ux{0) = '^^(O) for |a:| < K, for each x G \—K,K] and T > 0, for lower 
boundary values we have 

0 < Ux{t) — ^Ux{t) < P {Js [—Ki K\ for some s G [0, i]). (11) 

and the same is true for upper boundary values if we exchange Ux{t) and ^Ux{t). In particular, for 
fixed L, T, 


^Ux{t) —>■ Ux{t) as K ^ oo uniformly for x G [—L,L] and t G [0, T], (12) 

Proof. We consider lower boundary values but the proof is analogous in the other case. That Ux(t) — 
^Ux{t) > 0 follows from the discussion just above about monotonicity property of the limiting dual. 
The other inequality in (jlip follows from the inclusions 

{Js C [-K, K] for all s G [0,t]} H {Ct ~ u} C {^Ct = Ct} H [Ct ~ u} C {^Ct ~ u}. 

For the second statement m, notice that from |/o| = 1 the cardinality of the set 

{(x,rc) G H : (Xjin) G Is for some s G [0,t]} 

is dominated by a branching process Zt with Zq = 1 in which each particle gives birth to a pair of 
particles at rate a + 6 . Now, let Co = {(^^i ii^)} where x G {—L, L], and let t G [0, T]. Then, since the 
diameter of the set Js can only increase by M at each birth event, 

P(J, C for all s G [0,t]) > P {Zt < {K - L)/M) > P {Zt < {K - L)/M) 

= 1- P{Zt > {K - L)/M) > I - M {K - L)-^ E {Zt) 

= 1 - M (X - L)-^ ^ 0 as iiT^oo. 

This shows statement ()12p and completes the proof of the lemma. □ 

Lemma [TT] gives only rough bounds on the rate of convergence of ^Ux{t) as K ^ oo. With some 
extra work we get a better bound, that will be necessary in the proof of Theorem 01 Note the n 
below has nothing to do yet with the patch size N, though it will later on. 

Lemma 15 - For variable n > 0 and fixed ai,a 2 , let T := ailogn and L := a 2 iogn. Then, for 
each K > 0, there is 03 and no so that, for K = log n, 

\^Ux{t) — Ux{t)\ < n~^ for all n > no, \x\ < L and t G [0,T]. 

Proof. This follows from the first part of Lemma fTTl In order to find a good upper bound for the 
right-hand side of we proceed in two steps by hrst controlling the growth of the set Js in a 
short time interval and then at some regularly distributed times. 

Step 1 - Let Zt be as in the proof of Lemma [TT] and let A := a -|- 6 . Then, 

P(ZT>n"i^+^) = P{ZT>n^e^'^) = P {Zt > E {Zt)) < n"^ 
for all ,5 > 0. In particular, letting e := aiA -|- (I, it follows that 

P {B) > 1 — n~^ where B := {|/t| < rJ for all t G [0, T]}. 


(13) 
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In addition, letting Xt count the number of transitions occurring in the influence set up to time t, 
we have the following stochastic domination; 

P {{Xt+h — Xt>k} B) < P (Poisson {n^s) >k) for all f < f + s < T. 

Letting h = n~^ with <5 > e > /3, we deduce that 

P {{Xt+h — Xt > k} n B) < P (Poisson (n^h) > k) < P (Poisson (n^“^) > k) 

( 14) 

— < (n'^“'^)^ (1 — 

for any 0 < t < T — h. 

Step 2 - We now look at the process Jg at the times in 


5 := {0 < t < T : t = jh for some j £ N}. 

Let Iq = {(xjU))} for some |x| < L and Jg^x be the corresponding values Jg. Then, it follows from 
Lemma [TT] that, for all 7 > 0, there is c > 0 such that, for s <t, 

P {Js,x [-L - ct,L + ct]) < P {Jg^x [x - ct,x + ct]) 

< 2 P {Jg^x ^ {—oo,x + ct]) < 2e~'^^. 

Replacing t with T and letting Kq := L + cT , for every m > 0, 

P{Jt,x^[-Ko,Ko]) < n-^ for all t<T 

for some c > 0. In particular, for m > 5, 

P {Jt,x ^ [—KoiKq] for some t € S) < \S\n~^ < ailogn (15) 


for some c > 0 . 

Conclusion ~ Let k so that k{6 — e) > 6 and K := Kq + kM where M is the dispersal range. 
Using (fTT)l at all times t £ S together with the estimates (fTH|l and m, we find 

P{Jt,x C [-K,K] for all t £ [ 0 ,r]) 

> 1-P {B^) - P {{Jt^x \-K, K\ for some t £ [0, T]} n B) 

> P{B)-P {Jt,x t [-^ 0 , Kq] for some t£S)-\S\P {{Xt+h - Xt > k} n B) 

> 1 — n~^ — ai logn — IS”! (1 — 

> 1 — n~^ — ai log n (1 — 

Using Lemma dH we deduce that, for x £ [—L,L] and t £ [0,T], 

\^Ux{t)-Ux{t)\ < n-^ + ailogn(n -(™-^)(1 
To conclude, it suffices to take 5, then m, then k such that 

(3 > K and m — 6 > k and k{5 — e) — 6 > k. 


This completes the proof. □ 
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5. Coupling of the dual processes, collision estimates 

We now show how to couple the Wduals to the limiting dual. We begin by constructing the limiting 
dual a bit more explicitly. For integer k > 0, define the following random variables. 

• Let {an{k),ui{n, k),U 2 {n, k) : n > 0} := Poisson point processes with rate a, and a pair of 
independent Uniform (0,1) random variables attached to each Poisson event. 

• For 0 < \m\ < M, let {bn{k,m),ui{n, k,m),U 2 {n, k,m) : n > 0} := Poisson point process 
with rate 6/2M, and a pair of independent Uniform (0,1) random variables attached to each 
Poisson event. 

• Let {dn{k) : n > 0} := Poisson point process with rate 1. 

In case {x, w) € 1^- has personal label k, then 

• at time t = an{k), set It = It- U {{x, ui{n, k)), {x, U 2 {n, k))}, 

• at time t = bn{k, m), set It = It- U {{x + m, ui{n, k)), (x + m, U 2 {n, k))}, and 

• at time t = dn{k), set It = It- \ {x,w), 

with the labels of newly added points assigned as described earlier. 

We now construct a copy of the Wdual using the same random variables. To distinguish it from 
the limiting dual, we use the notation I^^, etc. We first divide the interval [0,1) into the N 
subintervals [(j — 1)/N,j/N). Then, to each point in I^, we add a fourth label, which we call the 
location label and is assigned as follows. 

In case (y, v) is a new point in I^ and location labels 1,2,... ,k are in use, 

• if there is {x, w) G with x = y such that both second coordinates v and w lie in the same 

subinterval of [0,1) then {y,v) is assigned the same location label as {x,w), 

• otherwise {y,v) is assigned location label k + 1. 

To construct the W-dual from the previous random variables, we then make a slight adjustment to 
the above transitions, namely, if {x, w) G I^ has location label k, then 

• at time t = an{k), set I^ = I^_ U {(x, ui{n, k)), (x, U 2 (n, k))}, 

• at time t = bn(k, m), set I^ = U {{x + m, ui{n, k)), {x + m, U 2 {n, A;))}, and 

• at time t = dn{k), all points in I^ with location label k are removed from I^. 

Also, when a new pair of points is added to I^, with ui{n, k) and U 2 {n, k) being the uniform random 
variables attached to the Poisson event occurring at time t, their parent labels (now a subset of the 
integers Z) consist of the personal labels of all points with location label k. 

Definition 16 (collision) - We say that a collision occurs if a newly created point in I^ is 
assigned the same location label as an existing point. 

Note that, in the absence of any collisions, each point’s location label is identical to its personal 
label, and so the A^-dual and the limiting dual coincide. In particular, it is useful to estimate 

:= time of the first collision starting from the pair (q = {(x, rci), (x, 

Since the time of the first collision starting from a single point is stochastically larger, the following 
estimate holds for such a time as well. 
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Lemma 17 - We have P <t)< + 1) 0 as — >■ oo. 

Proof. The idea is to first control , which is the number of points added to in the time 
interval s G [0, t] including the two initial points at time 0, and then show that, as long as this 
number of points is not too large, the probability of a collision is small. Noting that is dominated 
stochastically by a branching process Zt starting with two particles and in which each particle gives 
birth to two new particles independently at rate \ := a + b, for each 7 we have 

P(nf >e^*) < P{Zt>e^^) < E [Zt] < (16) 

In other respects, if a new pair of points is added at time t, 

P (collision at time t \ n^_ < n) < 1 — (1 — n/N){l — (n + 1)/A^) 

from which it follows that 

P (r^ < t I nf < n) < {1/N) Yli<n i (17) 

Combining (fT6]) - (fT7)l with n = e'^'^, we get 

P{t^ <t) < P (nf > eT"*) + P (r^ < t I nf < e^'*) 

< + (l/Ar)e2P. 

The lemma follows by taking 7 = (l/3t) In A^ in the previous inequality. □ 

6. Occupation density 

For the A^-dual, we define the function in the same way as (pt has been defined in ([^ but using 
the process C,i^ instead of Q- Then, it follows from the above that 

|0t(u) — $f^(u)| < P(r^<t) for any u G AT (18) 

since the limiting dual and the A^-dual coincide as long as there is no collision. In case has more 
than one point, we write that ~ u when the label of every point in is active at time t with 
respect to n. Now, define the functions 

(</>?(u))x := P(Ct ~ u I Co = {{(a:,uii),(x,ui2 )}}) 

where wi and W 2 are independent Uniform (0,1), and similarly using C^- Since the sets 

evolving from the two points (x, wi) and (x, W 2 ) do so independently in the limiting dual, it follows 
that (ptiii) = as the notation suggests. Moreover, 

|.i.?(u) - <1.7(u)| < P{T«<t) (19) 

which motivates our definition of earlier. Our next result is concerned with the occupation 
density of the N'-patch model, which we define as (t) := {1/N) ^t{x)- 

Theorem 18 - For all e > 0 and x G Z, we have 

P(|«"(«)-«{u"{0))),|>e) < 2t-"F{T"<(). 
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Proof. We first prove the result subject to the symmetry assumption 

ji)) = 1 for i = 1,..., A:) = P (t/o((x, o-(ii))) = 1 for i = 1,..., A:) (20) 

for every subset {ii, j2) ■ ■ ■ ,jk} C {1, •••, A^} and any permutation a of this set. Note that if the 
previous equation holds for the initial conhguration r]Q, then symmetry of the evolution rules implies 
that the same is true of r]t for all t > 0. Using the duality relationship for the patch model together 
with the previous symmetry assumption ()20l) with A: = 1, we get 

Ei^tix)) = Ex:.(x)=.^’(^t(x) = l) = iV^'f(u^(0)). (21) 

Using the duality relation and (f20]l with k = 2, we get 

E{Ct{xf) = Ex,y:7r(x)=7r(y)=x = ??t(y) = 1) = {u^ (0)) ■ (22) 

Combining (fT8]) - (fT9]l and (f^ - (l22]) . we deduce that 

Var(6(x)) = iV2($2,Af(uA^(0))-(cl.f(u^(0)))2) 

< - <(>f(u^(o)) + 

< 2N^P{T^<t). 

This, together with Chebyshev’s inequality, implies that 

F(K(()-W(u'"(0))),|>e) = P(|««(«)-£(a«(i))|>e) 

< Var (tt,^(t)) = {eN)~'^ Y&i {^t{x)) < 2e“^P(r^<t) 

showing the result under the symmetry assumption ()20p . To prove the result in the absence of 
symmetry, it suffices to symmetrize the distribution of the initial configuration r]Q without changing 
the distribution of ^o- If the distribution of i]q concentrates on coiffigurations with a hnite number 
of individuals, then the symmetrized initial distribution i]q is given by 

P{^o = ri') = n PiV0 = v) 

x-4'{x)j^0 

where corresponds to r]' and in the sum, ^ corresponds to r]. More general distributions can be 
symmetrized by first applying the above formula to the finite-dimensional distributions, and then 
taking a limit. □ 


7. Properties of the mean-field equations and proof of Theorems [TH2] 


Recall the mean-held equations ([3|): 


u 


I 

X 



for all X € Z. 


First we address existence of solutions. Clearly, the set /C = [0,1]^ is invariant. In addition, since 
the right-hand side of (l3|) is Lipschitz continuous with respect to the sup-norm ||u||oo = sup^, \ux\-, 
standard theory guarantees existence and uniqueness of solutions. 

From the stability analysis of the single equation @ and Corollary 1101 it follows that 
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• if r < 4 then for any u(0) € JC, Ux{t) —>■ 0 uniformly in x as f ^ oo, and 

• if r > 4 and inf^; tia;(0) > U- then Ux{t) uniformly in x as t ^ oo. 

Therefore the zero solution is stable when r < 4, and the positive equilibrium Ux = is stable in 
the sup-norm when r > 4. 

We now address expansion and retreat, defined in the introduction, which leads to more robust 
notions of stability. We begin with a lemma on wave fronts, where the set of wave front is 

W := {u ^ 1C : Ux > Uy for all x < y}. 

In what follows, the property stated in Corollary 1101 is called monotonicity. 

Lemma 19 - If u(0) G W then u{t) G W for all t > 0. 

Proof. For z € Z,we dehne the shift map r^u on 1C by {tzu)x = Ux-z- Then, writing the system ([3]) 
as u' = F{u), we easily check that F(r^(u)) = Tz{F{u)), from which it follows that 

v(t) = Tzu{t) solves dl]) whenever u{t) solves (|3|l. 

For a wave front u, one has r^u > u for all z > 0. If u(0) is a wave front, by monotonicity, it follows 
that, for z > 0, Tzu{t) > u(t), therefore Uy-z{t) > Uy{t) for all y. In particular, for a\\ x < y it 
suffices to take z = y — x to see that Ux{t) > Uy{t). □ 

Using Lemmas [12] and [T9l we can prove the hrst part of Theorem [2] 

Lemma 20 - Expansion and retreat are open conditions. 

Proof. By continuity of the solutions to ([3]) with respect to a and b, it is enough to show that 
if expansion occurs for a given value u then tii(ti) > u for some ti, and that if retreat occurs for 
given values ui,U 2 then < ui for some ti. 

Since W is invariant according to LemmafTOl and u l(x < 0) G W, if expansion occurs with rti(fo) = 
u then u{t) > ul{x < 1) and so r_iu(fo) > u(0). Iterating, for integer m > 1, we get r_mu(mto) > 
u(0), which implies that Uximto) > u for x < m. Now, 

• let u{t) solve dH) with n(0) = u, 

• let t 2 be such that u{t 2 ) = {u + 2u+)/3, 

• let c, 7 > 0 be as in Lemma fT^ 

• let ts be such that 2e~F^ = (tt+ — u)/3, 

• let m be large enough that m > cmax(f 2 ,t 3 )- 

Using Lemma [12] and monotonicity, we obtain ui(mfo) > (2u -|- u+)/3 > u. Since the proof for 
retreat is similar, it is omitted. □ 

With a bit more work we obtain Theorem [T] 

Proof of Theorem [T] - In the previous proof, we showed expansion for a given value u im¬ 
plies ui{ti) > u and by invariance of W, Ux{ti) > rti(ti) > u, for some ti and every x < 1. 
Let e := tii(fi) — ti > 0, and let u(0) := ul{[—L,L]). Using this and Lemma [T^ it follows that 
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for all L large enough, and then by symmetry, 

Ux{ti) > ri + 2e/3 for all —{L + l)<x<L + l 

Comparing to the increasing solution v(t) with v(ti) = (rt + 2e/3) 1(Z) and using again Lemma [T^ 
if L is large enough then additionally uo{s) > u + e/3 for ti < s < 2ti. Using monotonicity and 
iterating, we find that for integer m > 1, 

Ux{mti) > u + 2e/3 for all x € [—{L + m), L + m] and 

> u + e/3 > u for all s>(m + l)ti and x^[—m,m]. 

The conclusion of Theorem [1] is satisfied for xq = 1, 5 = e/Z and c = 1/fi. Since the proof for 
retreat is similar, it is omitted. □ 

To conclude this section, we show in the next two lemmas the two quantitative estimates that 
complete the proof of Theorem [2j 

Lemma 21 - Expansion occurs for M = 1 if a + b/2 > 4 and b > 8/9. 

Proof. Let u(0) = ul{x < 0), for u to be determined. Since W is invariant and u(0) G W, 
u(f) G VP for f > 0 which implies that > UQ{t) for t > 0. Therefore, uo{t) satisfies the 

differential inequality 

u'q > (a + 6/2) Uq (1 - uo) - uq. 

If a + 6/2 > 4 and u > 1/2 then liminft^oo uo{t) > u+, the upper equilibrium of ([1]) for f = a + 6/2. 
Since u+ > 1/2, let ti be such that uo{t) > 1/2 for t > ti, then for t > ti, ui{t) satisfies the 
differential inequality 

u'l > {auf + b/8){l — ui) — ui- 

Denoting the right-hand side P(ui), suppose that P{ui) > 0 for ui G [0,1/2]. Letting u{t) solve 
the equation u' = P{u) with u(0) =0, 

u_ < 1/2 < u{t 2 ) < u+ for some t 2 

where are the non-trivial equilibria of ([!]) for r = a + b. Letting u = u{t 2 ), u > 1 /2 and by 

comparison, ui{ti + 12 ) > u, which implies expansion. 

It remains to find conditions that guarantee P{u) > 0 for u G [0,1/2]. If a -|- b/2 = 4, 

P{u) = Q{u){u —1/2) where Q{u) := {—au^ + {a/2)u— b/4). 

For any u, provided a < 4b, 

Q{u) < Q(l/4) = a/8-6/4 < 0, 

which implies P{u) > 0 for u G [0,1/2). Moreover, P is non-decreasing with respect to a,b, 
and P(l/2) > 0 if a -|- 6/2 > 4. The intersection of the lines 

a-1-6/2 = 4 and a = 46 is at (a, 6) = (l/9)(32,8). 

By monotonicity we hnd that expansion occurs whenever a -|- 6/2 > 4 and 6 > 8/9. □ 
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Lemma 22 - Retreat occurs for M = lifa + 6<4 and b > 0. 

Proof. Given a, b with a + 6 = 4, we will produce a wave front solution u(t) with 

Ux{0) = u* > for X < 0 and U3;(0) = for x > L 

such that Ti(u(ti)) < u(0) for some ti. Iterating, for integer /c > 1, it follows that Tk{u{kti)) < u(0), 
which means Ux{kti) < u* for x > L — k. Letting k = L + 1 and using monotonicity, it follows that 
the solution with u(0) = u* l(x < 0) + tt* l(x > 0) has 

Ux{{L + 1) ti) < u* for all x > L — (L + 1) = —1, 

satisfying the conditions for retreat with tQ = (L -f 1) ti. 

Let f{u) = 4tt^ {1 — u) — u denote the right-hand side of (jl]) for r = a -|- 6 = 4. Then, 

• f{u) < 0 for u € [0,1] with /(O) = 0 and /(I) = —1, 

• f{u) has a local maximum at tt = 1/2 with /(1/2) = 0 and 

• f{u) has a local minimum at some u* G (0,1/2) with /(u*) < 0. 

Since M = 1, the mean-field equations are 

u'x = iaul +ib/2){ul_i+ul^^)){l-Ux)-Ux 
For e, L > 0 to be determined, let u* = 1/2 -|- e and 

• let Ua;(0) = u* for all x < 0, 

• let uo(0) = 1/2, 

• let Ux{0) be such that Uy_i + rty+i = for all y G [0,L), which gives 

Ux = -e{l + e)xf/‘^. 

where L is such that if ul-i,ul are decided by the above formula, then ul-i > u* > ul, 

• let Ux{0) = u* for x > L. 

Note that by choosing e > 0 small enough, we can make Ui_i(0) — u* as small as we like. For this 
initial data, we have the following properties: 

• u^(0) = f{ux) = fiu*) < 0 for all x < 0, 

• Uq(0 ) = /(1/2) = 0 by construction, 

• u^(0) = f{ux{0)) < 0 for X G (0, L) by construction, 

• < /(u*) + b for some 5 > 0 such that d —)■ 0 as e —>■ 0, 

• u'^(O) = f{u^) < 0 for all x > L. 

Taking e > 0 small enough that 5 < |/(u*)|/2, we have u/_^(0) < 0, and then redefining u* = 
1/2 -|- e/2, we still have u^(0) = f{u*) < 0 for x < 0. In addition, we have Uq( 0) < 0 since Uq 
is strictly decreasing with u_i. Since u^(0) < 0 for all x, and uniformly outside a finite set, we 
have u^(0) < —d < 0 for all x and some (probably different) 6 > 0. Since u(t) is strictly decreasing 
at time 0, for every s G [0, h] for some h > 0, one has u(/i) < u(0) and, by monotonicity, one easily 
shows that u(t -|- s) < u(t) for t,s > 0. Since is non-decreasing in Ux±i and is non-decreasing 
in Ux for Ux G [u*, 1/2], it follows that for all x G [0, L), we have u'^(t) < —b < 0 so long as Ux > u*, 
so there is ti so that Ux(ti) < u* for all x G [0,L) and in particular, ri(u(ti)) < u(0). □ 
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8. Block construction for expansion and proof of Theorem [3] 

In this section, we combine the behavior of the mean-field model with the estimates established in 
Section [6] to prove Theorem [3l To do so, we use the comparison technique described in [TJ section 4], 
also known as block construction, to couple the A^-patch model to an oriented site percolation model 
on the directed graph Ti with vertex set 

H := {{z, n) G Z X : z + n is even} 

and in which there is an oriented edge 

(z, n) {z', n') if and only if z' = z^l and n' = n + 1. 

See Durrett [5] for a definition and review of oriented percolation. In the vocabulary of [7], suppose 
we are given a /c-dependent percolation model with density at least, i.e., in which sites are open 
with probability at least, 1 — 7, and let Cq denote the cluster containing the origin, i.e., the set of 
sites that can be reached from the origin by an open path. We recall the key ingredients from the 
reference, whose statements have been specialized somewhat to fit our context. 

Theorem 23 ([7], Theorem 4.1) ~ We have P (|Co| < oo) —>■ 0 as 7 —> 0. 

Given an initial configuration Wq C {z : (z,0) €41}, let 

Wn := {z : there is an open path (y, 0) —>■ {z,n) for some y € Wq}. (23) 

The next result will help us to get a stationary distribution. 

Theorem 24 ([7], Theorem 4.2 with p = 1) - Let Wq = 2Z. Then, 

liminf„_>.oo P (0 € W 2 n) > 0 whenever 7 > 0 is small enough. 

Suppose we have a collection of configurations A such that the truth of the statement 7 € A depends 
only on the values {^(x) : 7r(x) € [—L,L]} and let 

Xn ■= {m € Z : [m, n) ^ H and rjnT G A} (24) 

where for y € Z is defined by ri{x,j) = r]{x — y,j), and where L and T give the appropriate 
space and time scales for the block construction. The sets in (I23j) - (j24p are usually referred to as 
the set of wet sites and the set of occupied sites at level n. The next theorem gives a sufficient 
condition for the set of occupied sites to dominates stochastically the set of wet sites. 

Theorem 25 m, Theorem 4.3) - Suppose {T]t)t>o Is a translation-invariant finite range process 
and for each po £ A there is an event G{po) such that the following comparison assumptions hold: 

• the event G{pq) is measurable with respect to the graphical representation of the process 
in the finite space-time box [—kL,kL] x [0, A:T] where k is the range of dependency of the 
percolation process introduced above, 

• we have the inclusion {po € A} n G{po) C {a'^^ pT £ A} f] pr £ A} and 

• the scale parameters L and T can be chosen such that P (G(r/o)) > 1 ~ 7- 
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Then, the process Xn dominates Wn provided Wq C Xq . 

With the previous three theorems and the results from the previous sections in hands, we are now 
ready to prove Theorem [3l 

Proof of Theorem [3]- Suppose a and b are such that expansion occurs in the mean-field equations 
and let u,L,xo,6,c > 0 with u- < u < u+ he as in the statement of Theorem [H Then, 

^^ 3 ;( 0 ) > u for all |x| < L implies that Ux{t) > u + 6 for all \x\ < ct — xq 

so that, letting T := (3L + xo)/c, we get 

cT — xq = 3L and Ux{T) > u + 6 for all \x\ < 2>L. (25) 

Now, for an interval [—K, K] C Z, let denote the restriction of the A^-patch model to this 

interval with lower boundary values, which is the process constructed from rjQ by ignoring the 
Poisson processes 

a„(x,y,z), 6 „(x,y,z), dn{x) such that 7r(x),7r(y) or 7r(z) ^ [-iiT,i^]. 

In the same way, define the restriction of the Wdual and let ^Ux{t) denote the restriction with 
lower boundary values as given by Definition [131 Also, let 

A := {r/ : ^{x) > uN for all |a;| < L}, 

where ^ is the mesoscopic configuration corresponding to rj. That is, A is the set of microscopic 
configurations with at least uN occupied sites at each location in the interval [—L,L]. One easily 
verifies that Lemma [T71 and Theorem 1181 hold for the restrictions defined above. Using in addition 
Lemma [m and ([25]l . we deduce that, for all 7 > 0, there is K large such that 

P ^rjT S A and ^rjT S A | S A) 

= P > uN for all |x| < 3L \ ^tjq G A) 

> P > N i^UxiT) - 5/2) for all |a;l < 3L | ^rjo € A) 

> 1-{3L + 1) x2{5/2)-^ P{t^ <T) > 1-7 

for all N large enough. This shows that there exists an event G{r]o) that satisfies the last two 
comparison assumptions in the statement of Theorem |25l Since the estimates in ()26p holds for the 
process where births outside the interval [—K, K] are suppressed, this event can also be chosen to 
be measurable with respect to the graphical representation in the corresponding space-time box, 
showing that the first comparison assumption is satisfied as well. In particular. Theorem 1251 can be 
applied which, together with Theorem [23l implies that there exists L such that 

lim 7 v->.oo P iCtix) = 0 for some t > 0 I 7?o € A) = 0 

By attractiveness, the same holds when starting with all patches in [—L, L] fully occupied. To also 
prove the existence of a nontrivial stationary distribution, we consider the process starting from 
the all occupied conhguration. Due again to attractiveness, the distribution of pt is stochastically 
decreasing in time and by the general result m Chapter III, Theorem 2.3], the distribution of pt 
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converges to a translation invariant, and on each patch also permutation invariant, stationary 
distribution v. Combining Theorems 124112^ we get 

v{{rj : r/(x) = 1}) > 0 for some x G 

but by translation and permutation invariance, the same is true for all x G In particular, v 
does not concentrate on the all-zero configuration and the proof is complete. □ 

9. Block construction for retreat and proof of Theorem |4] 

This section is devoted to the proof of Theorem 01 That is, we show that if the inner and outer birth 
rates a and b are such that retreat occurs in the mean-field equations, then for N large enough the 
distribution of the Ai-patch stochastic model converges weakly to the point mass on the all-zero 
configuration. The strategy is as follows. 

1. Wait until is uniformly small for all |x| < L, for some large enough L. 

2. Using a block construction, attempt to bootstrap from this configuration to percolation of 
low density regions in space and time. 

3. If the construction fails, then wait until a low density region appears and try again. 

4. If the construction succeeds, then show that there is a completely vacant zone that grows 
asymptotically linearly in time. 

Since for fixed L and N the number of sites x such that |vr(x)| < L is finite, for each h > 0, there is 
a positive probability that in the time interval [t, t + h] every site becomes vacant without any new 
births onto those sites occurring. Thus, although it may take a long time, for any stopping time t, 
the first time s such that ^t+six) = 0 for all |x| < L has an exponential tail and, in particular, is 
finite with probability one. 

The idea of retrying the construction until it succeeds appears in m and other references. If 
the construction succeeds with positive probability and fails within finite time when it does not 
succeed, then we need only to try a geometrically distributed number of times until it succeeds. 
Standard techniques (see for example the contour arguments in the appendix of m) guarantee that 
for a finitely dependent percolation model with density > 1 — 7, if 

|IUo| < 00 and P {Wn / 0 for all n) < 1 

then conditioned on the extinction event {Wn = 0 for some n}, the least such value of n has an 
exponential tail and in particular, is finite with probability one. 

Our first step is the following improvement of Theorem [TJ 

Lemma 26 - If retreat occurs, then there exist u,L,c > 0 with 0 < u < U- and, for all e > 0, 
there exists an xo(e) > 0 such that 

Ux(0) < u for all |x| < L implies that Ux(t) < e for all \x\ < ct — xo(e). 

Proof. Using Theorem [T] we have u,L,xi,6,c> 0 with 0 < rt < u_ so that 

'W3;(0) < u for all |a;| < L implies that Ux{t) < u — 5 for all |a:| < ct — xi. (27) 
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Let v{t) be the solution to ([4]) with u(0) = u — 6, and let co,7o be as in Lemma US Since v(t) 
converges to zero as t —)• oo, 

max (n(to), < e/2 for some to > 0. 

Let xo := cto + coto + xi. Then, according to (l?71) . for fixed t, 

|x| < ct — xq and |y — x| < coto imply that \y\ < ct + coto — xq = c{t — to) — xi 

imply that Uy{t — to) < u. 

In particular, letting v(s) be the solution to ([3]) with Ua;(s) = v{s — (t — to)) for all x, using Lemma[T2] 
and monotonicity, we conclude that 


Ux{t) < Vxit) + e/2 < e 


as desired. □ 

Starting from a low density region, we create a larger, completely vacant region in two steps in the 
following manner. First, we show that after some time the density becomes low and then remains 
low for a while. Then, we use the maintenance of that low density to show that the region becomes 
completely vacant. We begin with the first of these two steps. For ai,a 2 ,a 3 > 0 to be determined 
in a moment, we let 


To := a\ logA^ and L := a 2 logA^ and K := as logW 
Let u, T,xo(e),c > 0 be as in Lemma [26] and, in the terminology of Section [8l define 

A = {7 : ^(x) < uN for all |x| < L} 

where ^ is the mesoscopic conhguration corresponding to rj. We want to make the density small on 
some region, and then exploit this fact to show that the population goes to zero in that region. Since 
we need an upper bound, dehne the iL-restrictions using the upper boundary values ^Uy{t) = 1 
for \y\ > K, so that ^Ux{t) > Ux{t) for all x,t. For 0 < e < min (u, (a + 6)“^) to be chosen later, 
dehne 

Aq = {7 : ^t{x) < eN for all |x| < 4T} and Go{r]o) = {^rjt G Aq for all t € [To/2, Tq]} 
Then, we have the following estimate. 

Proposition 27 ~ For each e > 0 and small enough ai > 0, there are 02,03 > 0 and 6*1,71 > 0 
so that, for all N large enough, 

P{Go{vo) I ho€ A) > l-GiN-^P\ 

Proof. Let u,Lo,xo{e),c be as above and hx xq = xo(e/4) so that 

'Wx(O) < u for all |x| < Tq implies that Ux(t) < e/4 for all |x| < ct — xq. 
Monotonicity implies that if T > Lq, then 

'Wa;(0) < u for all |x| < L implies that Ux{t) < e/4 for all |x| < ct + T — (Tq + ^^o)- 
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Fixing 02 > 0 small enough that a\ > we have 

c(To/2) + L — (Lq + xq) > 4L when is large enough. 

Then, take 03 large enough so that, by Lemma [T5l 

^Ux{t) < Ux{t) + for every \x\ < AL and t G [0,To]. 

Now, combining Lemma [T71 and Theorem 1181 we get 

P{^^t{x)/N <^u^{t)+ e/A) > l-2(e/4)-2p(T^ < t) 

> 1 - 2(e/4)-2 (2e2(“+^)‘+ l)iV-V3 > 1 - 2 (e/4)-2 (2 lV2(a+&) ai ^-i/3_ 

for fixed |x| < K and t G [To/2, Tq]. Letting 71 := 1/3 — 2(a + b) oi, which is positive for ai small 
enough, and using the above observations, we deduce that 

P{^it{x)/N <^e/A) > l-6(e/4)-2iV-3'i = 1 - (QO/e^) . 

Suffering a factor of 8T + 1 = 802 log + 1 in the estimate, this bound is uniform over |x| < 4L. 
To make it also uniform in time, let h = (e/8)(a + 6 + 1)~^. Each transition in the process increases 
or decreases the number of occupied sites by at most one. Moreover, for each x, the total rate of 
Poisson point processes affecting sites in patch x is at most (a + 6 + 1) N". Thus, for each j, 

suPsG[ift,(i+i)h] I6(a;) -ijh{x)\ 

is dominated by X = Poisson (/i) where ^ = {a + h + P) Nh = eN/8. A standard large deviations 
bound gives P {X > 2/i) < e~^P. Summing over |x| < 4L and To/(2h) < j < Tq/Zi, 

P{^it{x)/N < e for all |x| < 4L and t G [ro/2,ro]) 

> l-(8a2 logAr +l)(ai/(2/i)logiV)(96/e2)lV-T'i 
- (802 log N + l)(ai/(2/i) log N) e-(“+''+i) 

for all N large enough, and the result follows. □ 

If we only wanted percolation of space-time boxes with low density then the good event Gq would 
suffice. In order to take low density to no density, however, we need a bit more work, and a couple 
more events. Let T := 2 logA^ and dehne 

Glim) = e Ao for all t G [ro/2,T + To/2]} 

Since e < u, we have Aq C A therefore, iterating at most Ajai -|- 1 times the proof of Proposition [27] 
by shifting things upwards in time by To/2 at each iteration, we find that 

P (Gi (r/o) I % e A) > 1 - (4/ai + 1) Ci 

Now, dehne the new collection of configurations 

Ai = {?7 : ^(x) = 0 for all |x| < 3L} 

and fix 77 G Aq. For t G [To/2,T], we define the following processes which are coupled by being 
constructed from the same graphical representation: 
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• let rj^ be ^rjt as defined above, started from rj at time To/2, 

• let 77/ be ^^r]t started from 77 at time To/2 with lower boundary values, 

• let 77/ = 77/ — 77/ > 0 be the occupied sites in 77/ that are vacant in 77/. 

By construction, 77/ is the set of occupied sites x with |vr(x)| < 4 T that would not be occupied 

without the help of an outside site, i.e., an occupied site y such that |vr(y)| > 4 L. Let 

G 2 ,o{v) = {^T+To /2 e Ai} and ^2,1(7?) = G Ai for all t € [To/2,T + To/2]} 

where ^l{x) := card {x : ' 7 r(x) = x and r/Kx) = 1} for 7 = 0 , 1 , and define 

G2(r?) = G 2 ,o(r 7 )nG 2 ,i( 7 ?). 

Now, we say that there is a type 2 active path (xo,to) —^2 (x, t) if there are 
xo,xi,... ,Xn_i = X G Tat and to < ti < ■ ■ ■ < tn = t 

such that the following two conditions hold: 

• for 7 = 0 , 1 ,... ,77 — 1 , we have rj^ixi) = 1 for all s G [tj,L+i] and 

• for 7 = 1,2,. ..,77 — 1 , there is a birth event at time ti where the offspring is sent to Xj and 

where Xj_i is one of the two parents’ locations. 

This active path is said to cross the space-time rectangle R = [^1,^2] x ['S'1,5'2] when 

Si < to < t < S2 and min ( 7 r(xo), 7 r(x)) < Ai and max ( 7 r(xo), 7 r(x)) > A2. 

Then, we define the event G3{rj) and the set B as 

Gsirj) = {no active path crosses [— 3 T, 3 L] x [To/2,T + To/2]} 

B = (x G Z : |x| < 3 L and > 0 for some t G [To/2, T -|- To/2]}. 

By definition, we have G2,i(77) = {B = 0}. Moreover, since the distance from the offspring to each 
of its parents is at most M, if the set B has a gap of size M, i.e., 

B n {x,x + 1 ,... ,x + M} = 0 for some (x,x-|- 1 ,...,x-|-M} C [— 3 T, 3 T] 

then G3{t]) holds. In particular, for N large enough, G2,i{'i]) C G3{r]). Finally, define 

Gim) = UnGAi Gi( 77 o) n 1^77^0/2 = r]}n ^2(77) n 03(77) 

where, after choosing a set of representatives, the union is really over a hnite set, as depends 
only on sites x such that |vr(x)j < K. One verifies that G{r]o) is measurable with respect to the 
graphical representation in the space-time rectangle 

[—K,K]x[ 0 ,JT] where J=l + ai/A. 


In addition, whenever 

T (Gi(77o) 1 770 G A) > 1 - ei and T (02(77) ] 77 gAi)> 1 - 


€2 
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after conditioning and reassembling and noting G 2 C G 3 , we find 

P (GiVo) I % G A) > 1 - ei - €2. 

Since > rj and using Gi, on G{rio) we have G A as required. 

To deduce the existence of vacant regions we pick apart Theorem [25] somewhat. Given rjnT, we 
define the event corresponding to G{r]o) but on the rectangle 

[—K + 2zL, K + 2zL\ x [nT, nT + JT] 

and similarly for Gi, G 2 and G 3 . Recall that Xn = {z G 2Z + n : a~‘^^^rjnT G A}. Letting Yq = Xq, 
we define inductively the sequence 

Yn+i = {^; G 2Z + n : (z ± 1, n) G Tn and holds} 

where the ± denotes or and not and. We note that the proof of Theorem [25] given in [7| proceeds 
by demonstrating the stated properties hold for the collection (Tn)n> 0 ) and then trivially deducing 
the same properties for (Ai„)„>o, so we may as well use Yn- For n > 0, define 

in ■■= infy„ and := supT^ 

where we use the convention inf( 0 ) = 00 and sup( 0 ) = — 00 . 

Lemma 28 - For all re > 0, we have 

^nT+To/ 2 {x) = D for all X e[{2£n-l)L,{2rn + l)L] (28) 

and for all n > 1 and f G [(re — 1) T + To/2, nT + Tq/ 2\, we have 

^tix) = 0 for all X G [(2£„ + 5) L, (2r„ — 5) Fj. (29) 

Proof. The proof is by induction on re. By construction, we have Yq = {0}. Moreover, 

G(r/o) = G*'°’°^(ryo) Folds and ii =—1 and ri = 1 whenever = {0} 

and one verifies from the definition of G that (1281) holds. Therefore, suppose that (I28p holds for 
some re. We show that if zi < Z 2 ^ Yn and and hold then 

iln+i)T+To/ 2 {x) = 0 for all X G [{2zi - 3) L, (2^2 + 3) L] (30) 

and for all t G [nT + To/2, (re + l)r + To/2], 

^t(x) = 0 for all X G [{2zi + 3) L, (2z2 — 3) Lj. (31) 

Since this is true for one of the four combinations zi = £n+i ± 1, ^2 = ^n+i ± 1, the result follows. 

Now, by assumption and since in Y zi < Z 2 < Xn, we have 

^„T+To/ 2 (a;) = 0 for all x G [{2zi - 1) L, {2z2 + 1) Tj- 

Also, by G^ 2 o’^^ and G^q^\ we have 

^(n+i)T+ro/ 2 (^) = 0 for all x such that |x — 2ziL\ < 3L for z = 1,2. 
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If zi > Z 2 — 3 then (IHOl) holds, and (I^TT) is vacuous, so we are done. Otherwise, if we had 

Ct{x) > 0 for some {x,t) G {{2zi + 3) L, { 2 z 2 — 3) L) x [nT + To/2, (n + 1)T + To/2] 

there would be an active path crossing [(2zi — 3) L, {2zi + 3) L] for either i = 1 or i = 2, which is 
ruled out by the event □ 

By Lemma [28l to conclude that percolation implies weak convergence of the distribution of to 
the all-zero conhguration, it suffices to show that if percolation occurs then in —>■ —oo and ^ oo. 
Since nothing is lost by increasing to the value uN on [—T, L] and setting it equal to one else¬ 
where, we may assume is symmetric about reflection across x = 0. In this case, the distributions 
of the left and right edges in and can be deduced from one another by symmetry, so it suffices 
to show that either in ^ —oo or rn ^ oo when percolation occurs. Suppose 

|^n| < x and |r„| < r for all n G N and some (random) r < oo. 

Then, for each n, since Yn C [—r, r], the event 

En := does not hold for any 2 ; G [—r, r]} 

has probability P (En) > e(r) > 0. Since the percolation model is Ai-dependent for some k < 00 , 
for Hm = {k+l)m the events {En„/)m>o are independent so, almost surely, after a geometric number 
of attempts, En„^ occurs and so there is no percolation. 

Since, as noted above, G 2 ,i C G 3 for all N large enough, it remains to estimate P (G 2 ), which is 
done in the next lemma. 

Lemma 29 - There are C, 7 > 0 so that P (G 2 ) > 1 — GN~'^. 

Proof. Since G 2 is the intersection of the two events G 2 ,o and G 2 ,i, it suffices to show the lower 
bound for each of these two events. 

The event ^2,0 - Let = X^|3,|<4i Q- Then {Sj, = 0} C G2,o and 

S/ :< Zt for all t G [ro/2,T -|-To/2] on the event Gi{r]o) 

where Zt is a branching process with .^ro /2 = '^Tb /2 — T ™ which each particle dies at 

rate one and produces single offspring at rate (a -|- 6) e < 1. Therefore, 

P{G 2 ,o) > P{St+To/ 2 ^^) - -P (-Z't+To/2 = 0) = I - P {ZT+To/ 2 >f) 

> 1 — E Zrpj,_j’^i 2 > 1 — (8T-|-1) eA^ exp(((a-|-6) e — 1) T) 

> 1 - (8a2 logiV-h l)eiV-i+2(“+'>)^ 

since T = 02 logA^ and T = 21ogA^. In particular, for e > 0 small so that 2{apb)e < 1, the 
probability of the event ^ 2,0 goes to one exponentially with N, as desired. 

The event ^ 2,1 - We start with the following two observations: 

• Outside sites, i.e., sites y with |7r(y)| > 4T, lead to occupied sites inside [—4T,4L] at combined 
rate at most 2M[a -|- h)N, and at locations x with 4L — M < |7r(x)| < 4T. 
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• Each of these, in turn, leads to occupied sites inside [—4L, 4L] at rate at most A = (a + 5) e 
and dies at rate one, and each offspring is at most M patches distant from its parents. 

We call the sites that become occupied in the first item first generation, and the ones in the second 
item the descendants. In order for G 2 ,i to occur it suffices that the subtree of descendants of each 
first generation site have depth at most L/M — 1. To estimate the probability of this event, we thus 
look at the number of hrst generation sites and the depth of a typical subtree. 

Size of the first generation - The size X of the first generation satisfies 

P{X > AM {a + b) NT) < P (Poisson (2M(a + b) NT) > AM {a + b) NT) 

< exp(-M(a + 6)iVr/2). 

Depth of a subtree - Letting fj, := X (1 + A)“^ where A = (a + 6) e, the descendant subtree of each 
first generation site is dominated by a Galton-Watson tree with geometric offspring distribution = 
(1 — /i) which has mean m = p(l — For such a tree, the corresponding discrete-time 

branching process Zj has expected value EZj = therefore 

P{Zl,m>^) = P(^l/m>1) < E{Zl,m) = 

The probability that a subtree has depth at least L/M is thus 
< exp (-M(a + h) NT/2) + AM {a + b) NT 

< iV-(“+ft)A^^ + 8M(a-F6) ln(7V)iVl-(“2/M)|lnm| < jy-l 

for all N large and e > 0 small, since lim£_).o m = lim£_).o (a + 6) e = 0. □ 

10. Block construction for spread and proof of Theorem [5] 

This section is devoted to the proof of Theorem [5l which follows the same three-step process as the 
proofs of the previous two theorems. More precisely, 

1. we hrst study the mean-held equations ([3]) starting with a single fully occupied patch and all 
the other patches empty, 

2. we then use the convergence in distribution of the Wdual to the limiting dual in order to 
show that, at least in a bounded space-time box, the stochastic process behaves almost like 
its deterministic counterpart when N is large, 

3. we hnally use a block construction to deduce that the probability of survival of the stochastic 
process starting with a single fully occupied patch approaches one for N large. 

Motivated by the monotonicity result stated in Corollary fTOl we hrst study the system where two 
adjacent patches can interact while all the other patches remain empty at all times. This system is 
used to understand a single time step in the block construction. That is, we let 

u = Ux and v = Uy where x, y G Z with |x — y| = 1 

and study the following system of coupled differential equations: 

u' = F{u,v) := {au"^ + {b/2) v‘^){l — u) — u 
v' = F{v,u) := [av"^ + {b/2)u‘^){l — v) — V. 


( 32 ) 
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Assuming that r = a + 6/2>4 and letting 

Li_ := 1/2 — (1/4 — 1/r)^'^^ and := 1/2 + (1/4 — 1/r)^/^ 
some basic algebra implies that 

F{u,u) = —ru{u — u-){u — u+). (33) 

To prove spread of the mean-held model, we proceed in two steps. 

Step 1 - Starting with patch x fully occupied and patch y empty, the density in both patches 
converges to the locally stable hxed point 

Step 2 ~ Assuming that the density in patch x is at least at all times, the density in patch y 
converges to as well regardless of its initial value. 

These two steps are established in the next two lemmas, respectively. 

Lemma 30 - Let r > 4 and b > 2. Then, 

u{0) = 1 and u(0) = 0 imply that limi_>.oo u{t) = limi_>.oo v{t) = u+. 

Proof. To begin with, we observe that, whenever u + v = 1 and b > 2, 

{u + vY = {au^ + {b/2)v'^) V — u T {av'^ + {b/2)u^)u — V 

= {a — 3b/2) uv + b/2 — 1 > min ((l/8)(2a-|-6) — 1, 6/2 — 1) (34) 

= min (r/4 — 1,6/2 — 1) > 0. 

In addition, using (p3]l and that v !->■ F{u, v) is nondecreasing, we get 


u' 

II 

> 

F{u, u) 

> 

0 

when 

U- < u < min (u, tt+) 


u' 

II 

< 

F{u, u) 

< 

0 

when 

U > V > u+ 

(35) 

u' 

II 

< 

F{u, u) 

< 

0 

when 

u + V > 1 and u > and v < U-. 



These three inequalities and their analogs obtained by switching the roles of u and v are summarized 
in the picture of Figure El Combining (|34p and the last statement in (|35l) implies that, under the 
assumptions of the lemma, there exists a time t finite such that 

{u + v){t) > 1 and min (tt(t), u(t)) > tt_. 

This, together with (|33p and the first two statements in (|35l) . implies the lemma. □ 

Lemma 31 - Let r > 4 and 6 > 2. Then, there are e > 0 and T < oo such that 

u{t) > u+ — e for all t <T implies that limi_>.oo u{t) = limt_>.oo v{t) = 

Proof. We prove the result when u{t) = u^ — e for all t < T and o < 8. In view of the monotonicity 
of the system, this will imply that the lemma holds in the general case. Under these specific 
assumptions, the second equation in (p^ becomes 

v' = {av'^ + {b/2){u+— — v) — V 

= (an^-|-(6/2) tt^)(l — n) — u-I-(6/2)(e^ — 2eu+)(l — n) (36) 

> (tt+ — v) Q{v) — be 
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fB U+ 




Figure 3. Summary of the inequalities in (I35II showing the sign of the derivatives on the left and schematic picture 
of the two steps of the block construction used to prove Theorem [5] on the right. 


where Q{X) := aX"^ + a {u+ — 1) X + {b/2) For all X G (0,1) and a < 8, 

Q{X) > Q{-a{u+-l)/2a) = Q((l/2)u_) 

= (a/4) + (a/2)(tt+— 1) u_ + (6/2) (37) 

= - {a/4)u‘i + {b/2)u+ > 6/4-a/16 > 6/4-1/2 > 0. 

Combining (l36]) - (|37|l . we deduce that, for all e > 0 small, 

V < 1/2 implies that v' = (tt+ — l/2)(6/4 — 1/2) — be > 0. 

In particular, for all e > 0 small, there exists T < oo such that 

u{t) = — € for all t < T implies that u{T) >1/2 and v{T) > 1/2. 

This, together with the direction of the arrows in Figure El shows that 

limt^oo u{t) = limt^oo v{t) = u+ 

This completes the proof. □ 

With the previous two lemmas in hands, we can now fix the appropriate time scale for the block 
construction. Let e > 0 be small. Then, according to Lemma [30l there exists a time Ti < oo, fixed 
from now on, such that we have 

v{t) > — e/2 for all t >Ti when u(0) = 1. (38) 
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Also, by Lemma EH there exists T 2 < 00 , fixed from now on, such that 

v{t) > U-^- — e/2 for all t > T 2 when u{t) > — e for all t < 2 T 2 . (39) 

For these deterministic times Ti and T 2 , we have the following lemmas which can be seen as the 
analogs of the previous two lemmas but for the stochastic process. 

Lemma 32 - Let r = a + b/2 > 4 and b > 2. Then, 

limw^oo Po (6(1) < (^^+ - for some t G (Ti, Ti + 2 T 2 )) = 0. 

Proof. Using Lemma [T71 Theorem [TSl and (E8]l . we get 

Po (6(1) < {u+ - e)iV for some t € (Ti,Ti + 2 T 2 )) 

= Po {u^{t) < — e for some t G (Ti, Ti + 2 T 2 )) 

< Po (|tt(^(t) — v{t)\ > e/2 for some t G (Ti,Ti + 2 T 2 )) 

< 2(e/2)-2p(r^ < Ti+ 2 T 2 ) < 8 e-^ {2 + 1) . 

This completes the proof. □ 

Lemma 33 - Let r = a + 6/2 > 4 and b > 2. Then, 

lim7v_^oo P (6(1) < iu+ — e)N for some t G {T2,3T2) \ 

6(0) > (u+ — e)N for all t G (0, 2 T 2 )) = 0. 

Proof. Following the proof of Lemma ES] but using (E^ instead of (ESl) . we easily prove that the 
conditional probability to be estimated is at most 

2(e/2)"2p(r^ < 3r2) < 8 e"^ (2 1) 

which goes to zero as —>■ 00 . □ 

Using the previous two lemmas, we can now prove the theorem. 

Proof of Theorem!^ - Declare site {z,n) G LT to be good whenever 

6 ( 2 ) > (u+- ei)A^ for all t G (Ti + nT 2 , Ti + (n + 2) T 2 ) 
and let Xn be the set of good sites at level n, i.e., 

Xn ■= {2 G Z : {z,n) G Lf is good}. 

Lemma E2] and obvious symmetry imply that 

Po ({-1,+1} iZlXo) 

= Po (min 2 =_i 4 ^t{z) < (u+ - e)iV for some t G (Ti,Pi + 2r2)) 

< 2 Po (6(1) 6 (ii+ — for some t G (Pi, Ti + 2r2)) ^0 as N ^ 00 


( 40 ) 
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while, according to Lemma [33] and again symmetry, 

P (z =h 1 ^ Xn+l I ^ G Xn) 

= P {^t{z ± 1) < {u+ - e)N for some t G (Ti + (n + 1) r2, Ti + (n + 3) T2) \ 

Uz) > iu+ - e)N for all t G (Ti + nPs, Ti + (n + 2) T 2 )) (41) 

= P (6(1) < («+ - e)N for some t G {T 2 , 3 T 2 ) | 

6(0) > (ti+ — e)A^ for all t G (0, 2 T 2 )) —0 as iV — 00 . 

Letting Wn be the set of wet sites at level n in an oriented site percolation process in which sites 
are open with probability 1 — 7 , the limit in (j41l) and Theorem 1251 imply that, for all 7 > 0, there 
exist N large and a coupling of the two processes such that 

Wn C Xn for all n G whenever Wq C Xq. 

This, together with the limit in (|40p . implies that 

Po (6 = 0 for some t > 0) 

< P {Xn = 0 for some n \ Xq = {-1, +1}) + Pq ({-1, +1} ^ Xq) 

< P {Wn = 0 for some n | Wo = {-1, +1}) + Pq ({-1, +1} Wq) 

which, according to Theorem 1231 goes to zero as N ^ 00 . □ 

11. Long range dispersal and proof of Theorem [6| 

This last section is devoted to the proof of Theorem [ 6 ] which shows that, for any choice of the 
inner and outer birth rates a and b and patch size N, the probability that the process starting 
with a single occupied patch survives tends to zero as the dispersal range tends to infinity. This 
together with Theorem [5] suggests that dispersal promotes extinction of processes with sexual 
reproduction, whereas dispersal is known to promote survival of the basic contact process with no 
sexual reproduction. The proof relies on the following two key ingredients: 

1 . In the absence of migrations: 6 = 0 , the process starting with a single fully occupied patch goes 
extinct almost surely in a finite time. This directly follows from the fact that, in this case, the 
process converges to its unique absorbing state because the state space is finite. 

2. Calling a collision the event that two offspring produced at the source patch x are sent to 
the same target patch y x, in the presence of migrations: 6 7 ^ 0 , but in the absence of 
collisions, offspring sent outside the source patch cannot reproduce due to the birth mechanism. 
In particular, in view also of the previous point, the process dies out. 

It follows that, to find an upper bound for the survival probability, it suffices to find an upper 
bound for the probability of a collision since we have 

Po (6 7 ^ 0 foi^ t > 0) < Pq (collision). (42) 

In addition, since the probability of a collision is related to the number of individuals produced at 
patch X and sent outside the patch which, in turn, is related to the number of individuals at the 
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source patch in the absence of migrations, to find an upper bound for the probability of a collision, 
the first step is to hnd an upper bound for the time spent in state j defined as 

r, := J^^P{Xt=j\Xo = N)dt for j = 1 , 2 ,..., 

More precisely, we have the following upper bound. 


Lemma 34 -Letb = 0. Then, Ej=i, 2 ...,N 3^3 ^ 

Proof. To simplify some tedious calculations and find the upper bound in the statement of the 
lemma, we first observe that the number of individuals at patch x is dominated by the number of 
particles in a certain simple birth and death process truncated at state N. More precisely, we note 
that the rate of transition j —)■ j + 1 is bounded by 

aj{j-l)N-\N -l)-\N -j) < afN-\N-j) < (a/4)j. 

In particular, standard coupling arguments imply that 

P{Xt>i\ Xo = N) < P{Zt>i\ Zo = N) for all i = l,2,...,N (43) 


where Zt is the continuous-time Markov chain with transitions 


j j + 1 at rate Pj := (a/4) j for j = 0,1,... , — 1 

j j — 1 at rate pj := j for j = 0,1 ,... , A^. 

Therefore, letting aj denote the amount of time the process Zt spends in state j, which can be seen 

as the analog of the amount of time Tj, inequality (1431) implies that 


Ylj=l,2,...,N 3 '^3 


= / r * . N P(X,=j\Xo = N) dt 

= X E ,.,2 . N Ej.,.A- P (A, = j I X„ = A) dt 

= /R,Ei.l,2. KP(X,>i\Xo = N)dt 

— / r + Ydi=i,2,...,N P [Zt > i \ Zq = N) dt = Yl,j=i,2,...,N 3 


(44) 


Now, to find an upper bound for the occupation times aj, we first let vj denote the expected number 
of visits of the process Zt in state j, which gives the recursive relationship 


Vj = Hj+i {/3j+i + fij+i) ^ Vj+i + /3j-i {/3j-i + fij-i) 1 Vj-i 
= 4:{a + A)~^Vj+i + a{a + 4:)~^Vj-i 


(45) 


for j = l,2,...,A^—1, with boundary conditions 


Vo = TiWi+Ti) ^^^1 and vn = I + /3n-i {Pn-i + Tn-i) ^ vn-i- (46) 


Note that the extra one in the expression of vn comes from the fact that Zq = N. We observe that 
the recursive relationship (I45p can be re-written as 

Vj = {I + a/4:)vj-i - {a/4)vj-2 for j = 2,3,...,N 
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which has characteristic polynomial 

(l + o/4)X + o/4 = {X -1){X - a/4). 

Using in addition that no = 1 since state zero is absorbing, we get 

ni = l + (o/4) and Vj = C 22 + 023 {a/4y = 1 + (o/4) H-h (a/4)-^ (47) 

from which it follows that 

Gj = = j-i(l + a/4)-i Ei=o,i,...j 

for j = 1, 2,... , — 1. Using ([MD-dU]), we also deduce 

VN-l) 

< ^ (1 + (o/4)(l + a/4) ^ 

— A^ ^ (1 + (a/4) (®/4)*) = X ^ 

Finally, combining (11^ with (I48|) - (|49]) . we get 

Si=l,2,...,Af U — Ylj=l,2,...,N 3 — Ylj=l,2,...,N Yli=0,l,...,j 

which completes the proof. □ 

Using Lemma 15^ we can now deduce upper bounds for the probability of a collision which, together 
with the inequality in (|42p . also give the theorem. 

Proof of Theorem - Let X denote the number of individuals born in patch x and then sent 
outside the patch before the patch goes extinct. The idea is to use 

Po (collision) = P (collision | X < M^S) (x < M^S) 

+ P (collision | X > MV3) (x > M^S) ( 50 ) 

< P (collision | X < (X > M^S). 

To estimate the first term in (j5Up . we observe that, since offspring sent outside the patch land on 
a patch chosen uniformly at random from a set of 2M patches, we have 

P (collision I a: < MV3) < 1 _ (1 _ j/2M) 

1 /■) (51) 

< 1 - (1 - MV3/2M)^ ^ < 1 - exp(-M2/3/2M) < (1/2)M-V3 

for all M large. For the second term, we first use Lemma (Ml to get 

= bN Ej=2,3...,N 3 (J - 1) N-^ {N - 1)-1 p 

— b Ylj=2,3,...,N 3 Tj — b J2j=l,2,...,N 


(48) 


(49) 


Px(^) 
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Looking at the different values of a, we deduce that Ex{X) is bounded by 

6 (1 — a/4)“^ 2 Af (1 “ when a <4 

^ Ej=i, 2 ,..,Ar (i + 1) ^ (V2)(^ + 2)^ when a = 4 
6 (a/4 — 1)“^ 2 AT («/4:y’'"^ < ^ (a/4 — 1 )“^ (a/4)^+^ when a > 4. 

This and Markov’s inequality Pq {X > M^/^) < Ex{X) give 

Po {X > MV3) < M-V3 5iv(i _ a/4)-i when a < 4 

< M-V3(5/2)(iV + 2)2 when a = 4 (52) 

< 6 (a /4 — 1)“2 (o/ 4)-^'^+2 when a = 4. 

Theorem [ 6 ] directly follows from (I50l) - (|52p also using (1421) . □ 
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